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We present a formalism for spatial averaging in cosmology applicable to general spacetimes and 
coordinates, and allowing the easy incorporation of a wide variety of matter sources. We apply 
this formalism to a Friedmann-Lemaitre-Robertson- Walker universe perturbed to second-order and 
present the corrections to the background in an unfixed gauge. We then present the corrections that 
arise in uniform curvature and conformal Newtonian gauges. 



I. INTRODUCTION 



The issue of the averaging, or fitting, problem in cosmology is an old one and addresses the question of whether a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metric is a reasonable approximation to the metric of the universe 
on large scales, and whether the evolution of the universe follows the usual FLRW dynamics. As the Einstein tensor is 
nonlinear, even should the metric of the universe be of FLRW form there is no guarantee that the same is true for the 
i— —\[ dynamics, and corrections therefore arise that act as an effective fluid. This effect is typically known as cosmological 
backreaction. The study of averaging in cosmology dates back at least to the 1960s [l[ and was revisited in subsequent 
decades (see for example @, y, [1, S S 0] for an inexhaustive selection of studies). Following the discovery of the "dark 
energy" in the late 1990s there was a great increase of study in this area (e.g. H, H, 0, El, E3 ) ■ Recent studies have 
i begun to quantify the effect in perturbed FLRW systems with varying results; see for example [H, [lj, [HI, [l6|, [ItJ 

and for further references see [l8|. However, many of these studies were limited to pure dust ( "Einstein-de Sitter") 
universes, and frequently to spatial slicings in s ync hronous coordinates. In previous work [l9[ (hereafter Paper I) and 
[HI (hereafter Paper II), as in for example [20l |2l|. a Newtonian slicing was employed instead. Paper II additionally 
introduced a general formalism applicable to systems with a variety of fluids for a wide range of redshifts. 

When considering the issue of backreaction in perturbation theory, the question of gauges has to be carefully 
, considered. The use of different gauges in backreaction introduces not least ambiguities in the definition of the spatial 
volume, which we consider particularly in flat gauge. Furthermore, it is not immediately obvious that there does not 
exist a gauge that is particularly suited to studies of backreaction - in which the correction terms take their simplest 
forms or, indeed, could be shown to vanish identically. (As an example, a flat gauge would remove the curvature 
| correction, while a Newtonian gauge would simplify the system by removing anisotropic terms.) 
• • ■ As yet it has been impossible to study this issue in its full generality, as to consider general gauges it is necessary 
to possess a formalism that incorporates a shift vector. More generally, a formalism that includes a shift vector can 
consider spacetimes in general coordinates, rather than those considering a particular symmetry. The central aim 
of this paper is to develop this formalism, and as a first example, apply it to perturbed Robertson- Walker universes 
in an unfixed gauge. 1 It is then a simple matter to choose two promising gauges - the uniform curvature, or flat, 
and the conformal Newtonian gauges - and reduce the general expressions to particular cases. We immediately find 
that, unsurprisingly, the backreaction cannot be removed through a judicious choice of gauge. Moreover, these two 
obvious choices remain relatively complicated, although uniform curvature gauge is arguably physically clearer when 
considering spatial volumes and Fourier space, while Newtonian gauge yields simpler forms. In Sjll]we present a brief 
overview of the 3+1 separation before generalising the familiar spatial averaging procedure to the most general case. 
We apply this general formalism to a linearly-perturbed FLRW universe in an unfixed gauge in mil Al and present the 
backreaction quantities in full detail before considering the specialisation to flat and conformal Newtonian gauges in 
ijlllBl We end with a brief discussion. 
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1 The issue of gauges in backreaction was considered in a study with somewhat different aims in |22l . while a formalism dealing with 
general coordinates was presented in |23t while the first version of this work was in final preparation. 
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II. GENERALISED AVERAGING IN COSMOLOGY 



The standard spatial averaging procedure in cosmology is based on a 3+1 split, where spacetime is foliated with 
purely spatial hypersurfaces upon which an averaging domain is defined. Assuming global hyperbolicity and a metric 
signature ( — h ++), a spacetime can be foliated with a family of 3-surfaces described by coordinates x^ = (t, x a ) and 
with a normal 4-vector 

n" = -(1,-/3*), n"n M = -l. (1) 
a 

Here the threading between corresponding events on neighbouring surfaces is defined by the lapse function a and the 
shift vector /3\ The projection tensor onto the 3-surfaces is 



and the induced 3-metric is 



= 9fj,u + n^n v (2) 
hj = g^Kh'j. (3) 



The line element is then 

ds 2 = {-a 2 + f3, t f3 l )dt 2 + 2[3 i dtdx i + h i3 dx l dx j , (4) 

where j3, — hijft . 

The embedding of the hypersurfaces is described by the extrinsic curvature, which is given by the Lie derivative of 
the metric along the normal, 

Kij — 2^ n ^ l, ij ■ ip) 

This can be rewritten as 

d 

aC n hij = T^htj - 2X>(i/3j) = -2aK tj (6) 

where T>i is the covariant derivative on the 3-surface and brackets denote symmetrisation on the enclosed indices. 

Using £>( ), and sffi respectively to denote a fluid's energy density, current density and stress tensor on the 
3-surface of a fluid species a, the 3+1 decomposition of the stress-energy tensor is 

T$ = Q{a)n»n v + 2n (M J^ + S$ . (7) 

Here 

Q (a) = n»n»T$, 4 a) = -n»h»T$ = h»4 a \ = h^T^ (8) 

and by construction n^J^ = n^sj^J — 0. We retain pt a \ and p/ a -\ to denote the fluid rest-frame energy density and 
isotropic pressure respectively. Projecting the Einstein equations onto the hypersurface and along its normal results 
in the Hamilton constraint equation 

K + K 2 - K)K{ = 16ttG ^ g {a) + 2A, (9) 

a 

where 1Z is the Ricci scalar on the surface and K = JQ, and the evolution equation 
d 



a dl 



Kn = a 2 f-Rij - 2K?K aj + KK %3 - teG^S® + AtirGhy £ (#(«) - g {a) ) - Kh l3 j 

\ a a / 

-a {V{D 3 a ~ (3 a V a K l3 - 2K a{z V ]) (3 a ) . (10) 



It also results in a momentum constraint equation with which we are not here concerned. For more details on the 
3+1 split, we refer the reader to [24l l25ll2(j|] . 
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The coordinate freedoms now allow us to eliminate four degrees of freedom. In previous studies of cosmological 
backreaction it has been standard to choose a synchronous slicing by setting a — 1, (3 l — (e.g. [j| [l(| EH). An 
alternative has been to employ a "Newtonian" slicing by choosing (3 l = and diagonalising the scalar sector of hij 
(e.g. [H, QjJ Hl| ) . In the following we do not fix our slicing and the coordinate freedoms are left implicit, allowing 
us to study spacetimes in arbitrary coordinate systems. 

Now let 



Vd= I Vhd 3 * (11) 
Jv 



be the definition of the volume of a finite domain T> lying on the inhomogeneous 3-surface where h is the determinant 
of the 3-metric. Then the average of a scalar quantity A in this domain can be defined as 



(A) = -1- / AVhd 3 *. (12) 
Vv Jv 



This averaging procedure cannot be sensibly applied to tensorial quantities since it is not invariant under a change of 
coordinates but retains its meaning for a scalar quantity. It has been used by most of the studies into cosmological 
averaging (e.g. 0, [TT|, 03); f° r a brief overview of alternatives see [13, HH. The averages are taken across spacelike 
hypersurfaces and therefore do not immediately relate to observation; for an approach at connecting a 3+1 averaging 
formalism with cosmology, see for example [29j ■ We can define an average "Hubble rate" from the expansion of the 
volume by 

3^ = (13 ) 
av Vv 

The derivative of the volume with respect to time is 

I = k L G^) ^ - k L { ~ aK + vlpi) ^ (14) 

where we have used that h = h^h^. The Hubble rate defined by the volume expansion is then 

3— = -<atf + X> i /3 i ) (15) 

with T> 1 — h^Vj. 2 Integrating this yields the average scale factor av- Combining the average Hubble rate with the 
averaged Hamiltonian constraint yields 

* a 

which can be interpreted as a Friedmann equation for ax> in the domain, where 

Kv = (a 2 K), Ql = (a 2 (K 2 -K i j Kl))-^(aK-V i f3 i ) 2 (17) 

are corrections from the averaged spatial curvature and the kinematical backreaction. The latter can be separated 
into the usual Qv (e.g. [1, HI, |T(| and Paper II) which is related to the variance of the extrinsic curvature, and 
which is dependent on the divergence of the shift, 

Qv = (a 2 (K 2 -K)Kl))~^{aK) 2 1 0% = | (D%) (2 {aK) - , Q£ = (18) 



2 The definition of the Hubble rate is not unique and it is common to instead define it from the expansion scalar of a fluid. The relation 
between the approaches was recently calculated in [23j . We retain our present definition to avoid connecting our slicing directly with 
fluid quantities while leaving free the possibility of doing so later. This simplifies the presentation at the cost of physical clarity. 
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Spatial averaging does not commute with differentiation with respect to time. The commutator between these 
operations applied to a scalar quantity A is given by 



(19) 



which can be written 

(A) = {AY + (A(aK- V l (3i)) - (A) (aK - V l p t ) . (20) 
The trace of the evolution equation for the extrinsic curvature is 



aK = a 2 + K 2 - \2ttG E g {a) + 4nG E S (a) - 3A j - aD^.a + af3 a V a K. 

\ a a / 



(21) 



Since 



(ak) = {aK)' + (a 2 K 2 ) - {aK) 2 + {aK) (V 1 ^) - (aKV 1 ft) - {aK) (22) 



and 

2 

a-D „ / a-D ' 



W = _ 3 ^ + 3 p +(&%)■, (23) 
we have after some manipulation and the use of (p~6|) that the average of pTjl can be written 



a-D 47rG 



E (« 2 tew + %))) + J (« 2 ) A + \ (Ql + + fe) . (24) 



a-p 3 
The dynamical backreaction [H, 0, [13] is 

= {aV l V ia ) - {aK) , (25) 
and the component arising from the variation of the shift is 

0% = + (X>*A) 2 - 2atfX> i /3 i - a^V.K) . (26) 

Equation (|24|) can be interpreted as a "Raychaudhuri equation" for the scale factor a-p. 

The energy density and isotropic pressure are defined on the 3-surface in l[7|). However, it is usual in cosmology to 
separate the stress-energy tensor with respect to the fluid's 4- velocity u^, 

T^ a) — n, My(°)«W j- „, , p( a ) _l 2o (a) 7/ ( ' a ^ 4- 7r (a ) -: — s n, \ — 7/^ t(«) 3 n , , — pernio.) (07) 

where 

= + "S„«f.) (28) 

projects tensors onto the 3-surface orthogonal to the 4- velocity, the 4- velocity is normalised to u^u^ — — 1, and 
pi a \ and p( a ) are the fluid's rest-frame energy density and isotropic pressure respectively. Let 

4°) = ^ ((> - ift) T&) = ((W - u U u U ) TV) (29) 

correct between the 3-surface and the fluid rest-frames. The effective, volume-averaged Friedmann and Raychaudhuri 
equations can therefore be written 



old \ 2 8irG 



a-D J 3 

&d 47rG 
od 3 



E (* 2 p(«)) + 1 <« 2 > A - 1 (**> + Qv - 6 E4 a) ) , 

a \ a / 

E (« 2 fa.) + Ma))) + \ (« 2 ) ^+\(VV + (% + Ql 3E4 Q) ) 
a \ a / 



(30) 
(31) 
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We therefore now have definitions of effective, averaged Friedmann and Raychaudhuri equations valid for any spacetime 
that allows a 3+1 split, with arbitrary lapse and shift, and an arbitrary number of fluids. Equations ([3H1 l3"Tj) are the 
main results of this work. 3 

It is tempting to immediately define an effective energy density, pressure and equation of state of the "backreaction 
fluid" and such can certainly be done. Simply interpreting ^ a (a 2 p( a )) and ^2 a (a 2 p^ a )) as the average density and 
pressure in the domain, the effective energy density and pressure of the backreaction can be defined as 

8nG v {a) 1, T 8nG 1 6 ffi + K D - AV V 4/3g - 3gg 

a 

with the effective equation of state 

Pcff 1 n v - 47> D - 4/3g - 3Ql + 6 J2 a 4 a) 

In some situations we should, however, take care with the interpretation of the results. A true Friedmann model is 
specified by a matter density, a cosmological constant and a curvature term. Depending on the time coordinate a 
conformal factor can also appear. Likewise, the corresponding Raychaudhuri equation is specified by a matter density, 
a matter pressure and a cosmological constant. Any or all of these terms can emerge from averaging the underlying 
structure: the effective cosmological constant would be a constant term that appears in the effective Friedmann 
equation, and the effective curvature term would be a term in the Friedmann equation proportional to 1/afj- The 
effective fluid might then sensibly be chosen so as to isolate these. Choosing the prescription l|32p would instead 
include the effective curvature in the effective fluid. 

When dealing with a true inhomogeneous model, this distinction is a matter of convention; the behaviour of 
the averaged scale factor remains the same. However our approach, as in Papers I and II, will involve specifying a 
background FLRW model and evaluating the corrections to this model from perturbations. The background Friedmann 
and Raychaudhuri equations will therefore in general include curvature terms and conformal factors. The averaged 
Friedmann equations, on the other hand, include these terms within TZx> and "Pp. For example, equation (3.9) of Paper 
II shows 1Zt> in a curved FLRW universe perturbed to second order. This includes the background term JC/a 2 . In this 
context - comparing the averaged equations to the assumed background equations - this term is not a correction. In 
other contexts this would not be an issue, but here it might be better to define the effective fluid in a model-dependent 
manner. 



III. BACKREACTION IN PERTURBATION THEORY 



A. Basic Results 



The perturbed Robertson- Walker line element with a scale factor a = a(rj) can be written as 

ds 2 = a 2 {r)) (-(1 + 2<j))drf + 2B l dr ] dx l + (5 y + 2C lj ) dx i dx i ) (34) 

where the perturbation <f> = ^2 n (J-/n\)(j)^ contains contributions from all perturbative orders, with corresponding 
definitions for Bi and Cy. In this work we only consider scalar perturbations, although the results of this section are 
equally valid for spacetimes containing vector and tensor perturbations. The spatial perturbation can be expanded 
into a trace and traceless component, 

Cy =-^5 ij +E ij , (35) 

where El = 0. In general we will only retain quantities to second-order in the perturbations and neglect all higher- 
orders. 4 To second order, the perturbations are 

= 0(i) + ^(2), B i = BW + ±B]?\ Ci^Glf + ^clf, (36) 



3 These equations can be compared with equations (36) and (37) of [23ll , 

4 This includes terms of the form (# 2 ) • In principle, these quantities could be large; however, it is clear that, within the confines of 
perturbation theory, these quantities are approximately of fourth-order in perturbations. 
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but except in products we never explicitly expand them. Indices on the perturbations are raised and lowered with 
the Kronecker delta. Derivatives with respect to the conformal time will be denoted with an overdot and H — a/ a is 
the conformal Hubble rate of the underlying model. Directly associating the spatial coordinates in (|34j) and (j4j and 
noting that f3 l — f3j gives 

h l2 = a 2 (Stj + 2Cij) , h l] =a- 2 (ft -2C lj + 4C* ifc C^) , fa = a 2 B„ f IV 2C ij B j . (37) 

Since a 2 — fa 1 fa = a 2 (I + 2</>), the lapse is then 

a 2 = a 2 (1 + 2(f) + B 2 ) (38) 

where B 2 = B l Bi. Paper II contains some discussion on the use of perturbation theory in backreaction, but some 
further comments are useful here. In essence, FLRW cosmology involves the use of an implicit "average" background 
upon which perturbations of arbitrary order are defined. When we apply a spatial averaging formalism to a perturbed 
FLRW cosmology, we are averaging across the perturbed manifold. Here, we are taking an average across a surface 
defined by perturbations that include second-order terms. This naturally introduces an apparent ambiguity between 
the Hubble rate a /a appearing in the usual background equations 

a\ 2 8nG 9 1 9 . a 4irG 9 . „ . 1 9 . ( a\ 2 

a 2 p+-a 2 A, - = —a 2 (p + 3p) + -a 2 A+ - 39 

3 a 6 6 \a J 

and in the averaged equations (|30l l3Tj) . However, the Hubble rate cl-d/cl-d is the Hubble rate defined by the time 
variation of an expanding volume living on the perturbed 3-surface, entirely distinct from the Hubble rate of the 
underlying model. In this paper, as in Paper I and Paper II, we will define the effects of "backreaction" to be the 
difference between these Hubble rates, and between the acceleration rates a-p/a-p and a /a. 
From the Ricci scalar on the three-surface, the curvature correction is 

Tip = ((1 + 20) {2d l &C tj - 2d l d t C) + AC l3 d l d j C + 4C lJ d k d k C l3 - 8C lj d k d l C{ + W l C jk d l C 3k ^ 

- (2 (0*CP*) (3^) ~ (d l C) (2^Gy - Z C) + 2 (di&t) {2d k C jk ~ djCf) ) . (40) 
The extrinsic curvature is 

aK) = --5} + \ (d l B 3 + d 3 B l ) - (d l C k + d 3 C lk - d k C}) B k + 2C ik C jk - C lk (d k B 3 + d 3 B k ) (41) 
and so the total kinematical backreaction is 

Ql = -4- (d'B, - B i {2&C l3 - diC) - 2("u),l>>, - \ (oV 

G 2 - Cpi + (d'B,) 2 - i (0*5,) (PBi) - l - (d l B 3 ) (BiB*) - 2Cd i B i + 2C)d l B^ . (42) 

The dynamical backreactions become 
a 



V v = 3 ( ~ ) + \ d * d ^ + ~ + C ) - -° lB i 

B,0> 0;W + {iriS,) (0,1V ) - (<9>) (di<f>) - 2C' 1 (),(),<:> - (<9>) (2d j C ij ~ d l C) ■ o (c <)' II, 
' WiB 1 + 600 + 2( •''('; , - 2C ij d i B j - B l (2d j C ij - d t C) ) (43) 



a 
a 



and 



13? = (d t B i + G-d i B i \- /2C ij d i B j +2C ij d i B j + l2^C ij d i B j + (d i B i 



- ( B* (2d j C ij ~ diC) + B l ( 2d j C ij - fyC) + 6^ {2d b C l3 - d l G) ) (44) 



B l (diC - di&B^ + 2Cd l B t -3^0% 
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The dynamical backreaction then includes a background term 3(d/a) 2 which appears in the background Einstein 
equations (|39p . To ensure that the "backreaction fluid" contains only correction terms, we include this in the averaged 
Raychaudhuri equation separately from Vt>- Since there is no gauge choice in which B a — C ab — 0, from (|40l l42l l43l 
I44[) it is already apparent that there exists no gauge in which all the correction terms vanish. 

Finally, we have the fluid corrections. For a perfect fluid with T$ = (p( a ) + P(a)) u ^ u * + PVaSd^vi the correction 
term is 

4 Q) = ^ (fa.) +P(a)) ((»M 0) ) 2 - « 2 )) = (iPM +P M ) ~ P l u[ a) ) 2 ~ a 2 )) • (45) 

The rest-frame energy density and pressure are expanded to second-order as 

P(a) = P{a) + S P(a) = P( )(l + \a)), P(a) = P(a) + S P(a)- (46) 

Here p^ = u>( a )/?( a ), and if we define the speed of sound to be c 2 , a s = dp( a )/dp(a) then the second-order expansion 
of the pressure perturbation is 

, dp( a ) l9 2 P(„), r . 2 _ / 2 s ( C s.(a))' r2 \ , An\ 

SP(a) = dp^ 6Pia) + = (5 W < a) " 2(l+t (a) )if ^ J ( 47 ) 

where 5r a \ includes contributions from all perturbative orders and we have used the background matter continuity 
equation. With the 4-velocity 



uf o) = a-Ml- 



(l - <t> + \ («?«) + 30 2 - 2B lV \ a) ) , v\ a ^ (48) 



where vf a) = <%^ a ) W W' J$ becomes 



(a) _ 8vrG 2 (- 



(P(a) +P(a)) ((v<i a) +B^) (y\ a) + S| o) )) . (49) 

Expanding (a 2 /5( a )) and (a 2 P( a ) ) with |J3BJ) and |@BJ>, we can finally express the Buchert equations in a perturbed 
FLRW universe in an unfixed gauge as 



dp \ 2 87rG 



« 2 E ?w + r 2 A - J + 2? - e e - 

a \ a / 

/ : S - E(^) + 3P(a))+^ 2 A+({ 



(50) 



d-p \ 4-7rG , 



| + # + Ql + 3C V - 3 £ Qt>> + 4 Q) ) ) 



(51) 



where now 



(a) _ 8?rG 2 



a-p(a) (*(«) + 2 ^ + 2^(a) + £ 2 + (1 + «> (o) ) (^ (a) + B() (y\ a) + B l ) ) , (52) 
•Sp ) = ^^a 2 P( a ) (3c 2 <5 (a) + 6w</> + 6c 2 (a) <M(a) + 3w (a) S 2 

+ (i + »«.,) («<" + Bf) („;„ + b U ) - 2(1 ( ft a ( '' ))g ^) (») 

correct the energy density and pressure respectively and £x> is a correction dependent on the cosmological constant, 

Cv = \a 2 K{2(j ) + B 2 ) . (54) 

The effective fluid giving the differences between the base and re-averaged FLRW model then has the effective energy 
density and pressure 

SrrG , 6^*^ + 6C V TZ V - Ql 

~~^~ a PcS — a i \ 00 ) 
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— « Pcff - 3 g (56) 

and an effective equation of state 

1 K V - 4V V - 4/3g - 3Ql + 12 £ a S { « ] - 18£ p 

U> ff = 7 — r . (£•') 

3 ^ + Qg-6E a 4 a) -6£ I , 

These can be compared with equations (33-35) of Paper II, generalised to a non-vanishing shift vector and conformal 
time. 



B. Gauge Choices 



The equations presented thus far are in an entirely unfixed gauge and can be adapted with relative ease to any 
perturbed FLRW universe. From the form of the correction terms in general, it is immediately apparent that we 
cannot find a gauge in which each contribution vanishes. The remaining question is whether we can find a gauge 
in which the corrections take on the simplest possible form. The most obvious choices are those gauges employed 
in standard cosmological perturbation theory for their simplicity: the uniform curvature, or flat, gauge; the uniform 
density gauges; the conformal Newtonian gauge; and the conformal synchronous gauge. Uniform curvature gauge 
would remove the curvature correction, a uniform density gauge should simplify the fluid corrections, the conformal 
Newtonian gauge should simplify the anisotropies significantly, and a synchronous gauge simplifies the dynamical 
backreaction. The gauge transformations and governing equations in FLRW universes perturbed up to second order 
can be found in [3l| and we will not repeat them here. 

The uniform curvature gauge is of immediate interest in backreaction. In this gauge the spatial surfaces align 
with those of the background FLRW model and the inhomogeneities are carried entirely by the choice of threading. 
As a result, by definition the curvature correction TZ-p vanishes in this gauge, although this is balanced against the 
introduction of the additional dynamical backreaction (3^ and shift terms in the kinematical backreaction. No other 
trivial choice of gauge entirely removes one of the backreaction terms when the system is written in conformal time; 
even for conformal synchronous gauge there is a non-zero contribution to Vt> arising from the time-derivative of a. 
A further benefit of the uniform curvature gauge is that for scalar perturbations the metric determinant reduces to 
\fh = a 3 , implying that the comoving volume of a spatial domain remains constant, which simplifies the averaging 
procedure significantly. 5 

The conformal Newtonian gauge is also easily motivated - while it doesn't remove any of the individual correction 
terms, the perturbation to the 3-metric is diagonalised and the shift is vanishing, which in particular considerably 
simplifies TZx> from the general case. It is also straightforward to find analytical solutions for simple systems in 
conformal Newtonian gauge, justifying its study (see for example [33l. l3~j| ) . 

Another interesting choice might be a uniform density gauge. Here the slicing is defined by choosing Sp = 0, and the 
threading by setting for example E = in (|35p . In this gauge, then, the velocity, shift and lapse are all non- vanishing. 
We can immediately see that the curvature correction will take the same form as that in conformal Newtonian gauge 
but that the kinematical and dynamical backreactions are significantly more complicated than in cither of the other 
two gauges. Furthermore, the aim of a uniform density gauge would be to simplify the fluid corrections 7r> and St>, 
and this it fails to do. Setting <5 = leaves the fluid corrections as unwieldy combinations of the lapse, shift and 
velocity. 

The final obvious choice is conformal synchronous gauge. We do not consider this gauge here, since the curvature 
correction lZz> is complicated (in real space) by the anisotropic components of the metric perturbation, and it has 
been well studied in the literature before (as for example in [l3l liH]). However, conformal synchronous gauge does 
have the advantage that the scale factor a d corresponds to a volume containing constant mass and therefore perhaps 
has a clearer physical meaning than in the other gauges. 6 

We will hence focus on the two choices that seem most likely to yield simple results: the uniform curvature and 
the conformal Newtonian gauges. In the current work it is not our aim to solve the resulting equations. Instead, we 
consider the form of the corrections in two particular, simple gauge choices. We will, however, briefly consider pure 
dust universes in both cases. 



In the terminology of [32l | and other works by these authors, a flat gauge is an example of a "volume-preserving" gauge. 
The authors are grateful to Dominik Schwarz for highlighting this point. 
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1. Uniform Curvature Gauge 

The uniform curvature gauge is of direct interest since the spatial hypersurfaces align with those of the FLRW 
background. The volume element on the hypersurfaces in flat gauge, and therefore the spatial averages across those 
surfaces, then coincide with the FLRW background. 

The metric in the uniform curvature gauge is found by setting C a b = in the above. Then 

ds 2 = a 2 {ri) (-(1 + 2(j>)drf + 2B i dr]dx i + <%cfa:W) (58) 
and so the spatial metric, shift and lapse are 

//,, </%. P 1 = B\ a 2 = a 2 (l + 2(j) + B 2 ). (59) 
This gauge being spatially flat, the curvature correction is trivially 

K v = 0. (60) 

The kinematical backreaction (l42l) becomes 



Ql = {{VB,) 2 - \ (d i Bj) (diB*) - \ (d 3 B l ) (BiB*) - ^-/bA . (61) 

Subtracting off the Hubble term, the dynamical backreaction (|431 Ef4")) is 

Vv = (d%4> + tf&djBi + (d l B 3 ) faBi) - (d^) - (j&Bi) - - (&Bi - 30 - W t B l + 6#) , (62) 

/?£ = (tfBi- (&Bi) 2 - tfdi&'B^ +3-(28 i B i - B'd^) . (63) 
In flat gauge the fluid velocity takes the form 



u, 



r 1 ( \ ... 1 / 



(«) = a 1 - + 2 + 2BtV ^ + W ) ' V ^ I (64) 



and so the fluid corrections become 



n-( a ) _ 



8ttG 



a 2 p(a) (*(«) + 20 + 2#(a) + B 2 + (1 + w (a) ) (« 4 (a) + Bi) (v^ + B^ ) (65) 



and 



$V } = ^ a2 P(a) (^ c l(a) S (a) + 6w (a) + 6c^ (o) #( a ) + 3w (a) B 2 ^ 

+ ((i + »<.,) («, w + *) (»:., + *) - Jl^c,) ■ < 66 > 

The correction arising from the cosmological constant is 

C v = \a 2 k(2^ + B 2 ). (67) 

We can now consider a pure dust (or "Einstein-de Sitter") universe, in which w — c 2 — A = 0. This case has the 
dual benefits of being one of the simplest cosmologies, and of being a reasonable approximation to the real universe 
for a range of redshifts from z w 1000 to z ~ 1. In this universe, the fluid corrections simplify slightly to 

T v = ^-a 2 p(S + 2^ + 2cl ) S + B 2 + (v l + B l )(v t + B i )), S v = ^^-a 2 p((vi + S f ) (V +-B*)) , £p = 0. (68) 
Working in Fourier space simplifies these expressions somewhat. Adopting the Fourier convention 

A k = y A(x)e lk x d 3 k, A(x) = y ^ke- 4k x ^, (69) 
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the spatial average of a quantity A(x) can be written as 

/H 3 k 
WUk)A(k)j^ (70) 

where W(x) is a window function defining the domain and W-p(x) = W(x)/Vx>. Let A(x) = /(x)g(x). If /(x) and 
<7(x) are both linear perturbations, we can take ensemble averages (denoted with an overbar) of both sides and use 
the definition of the power spectrum, 



2n 



2 



f{k)g*{k>) = —V{k)5{k-k') (71) 



to find 



1 r rjh 

(/(x) ff (x)) = ^3 J nk)f(k)g*(k) T (72) 

where we have used that a 3 W-p(k = 0) = 1. Writing the shift and velocity as 

Vi = ikiV, Bj = ikiB (73) 
therefore allows us to write the kinematical backreaction as 



In uniform curvature gauge, then, the entire kinematical backreaction is carried on averages of pure - i.e., non- 
quadratic - perturbations, and assuming these to vanish as in Papers I and II would unsurprisingly lead to an 
incomplete answer. This is the most important result of the calculation in this gauge. The other correction terms 
become 



V v = (J ( 3 ^( k ) - 4^fcB(k) - k 2 (j)(k)^ e~ ik - 



d 3 k 



(27T) 

' 1 ( - a (3B{k)B*{k) - 6<Kfc)</>*(fc)) - k<p(k)B*(k) - k 2 \<j>{k)\ 2 ^J V{k)^, (75) 



(2tt) 3 



Tv . (£) 2 (/( J(k ) + 2 *( k ))^ 

' " /" (20(fc) ( 5*(fc) + 2«(fc) J B*(fc) + 2|B(fc)| 2 + |w(A : )| 2 )7'(fc)*, (77) 



(2tt) 3 a 



5p = (~) J (\B(k)\ 2 + \v(k)\ 2 + 2v(k)B*(k))r(k)f. (78) 

While analytic solutions in this gauge can be presented, they lack the clarity of those performed in conformal New- 
tonian gauge, to which we now turn. 



2. Conformal Newtonian Gauge 

In Paper I and Paper II, along with for example [2l[, the Newtonian gauge was employed, as an alternative to the 
common use in studies of averaging of synchronous gauge. Newtonian gauge has certain advantages: the gravitational 
potential and spatial curvatures have clear physical interpretations and, moreover, remain small across all the scales 
studied. The Newtonian gauge also possesses a vanishing shift and a diagonal spatial perturbation, which simplifies the 
calculations somewhat. Here we use conformal Newtonian gauge. The major advantage in doing so is the simplicity 
with which solutions may be found - in Paper II the analytic solutions employed were derived in conformal Newtonian 
gauge. 
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Conformal Newtonian gauge is found by substituting <f> — > Bi — > and C a b — * —&8ab in the general expressions, 
and the line element is 

ds 2 = a 2 {r]) (-(1 + 2*)dry 2 + (1 - 2®)6 ij dx i dx j ) . (79) 
The kinematical backreaction becomes 

Ql = 6 (<i> 2 ) . (80) 

The dynamical backreaction is 

V v = (d i d l ^> + 3- (4> - $)\ + (2$<9 i a,* - (5% + <9 l $)) - 3 + 2-1"i< + 2-$$^ (81) 
and the curvature correction is 

K v = 4 (<9 4 <9 4 $) + (6 (<9 4 $) (9;$) + 8(* + 2$)<9 l <9 4 $) . (82) 

The fluid corrections are 

(a) _ 8nG 

and 



n -~"P(a)(5(a)+W + W6 {a) + (l+W {a) )vf a) ) (83) 



5 W = ^a 2 p (a) Ac 2 (a) <5 (a) + 6«; (o) * + 6< (o) <M (o) + (1 + ™ (q) > 2 q) - ^i^L- ^ g) \ , (84) 

and the contribution from the cosmological constant is simply 

£ C = \ 2 A(*). (85) 

These closely resemble those presented in Paper II, with only Vt> exhibiting minor differences from the change in time 
coordinate. 

In the particular case of a pure Einstein-de Sitter universe, the vanishing anisotropic stress implies that $ — ^ and 
that <i> = 0. The fluid corrections become 

Tv = ^a 2 p(5 + 2V + 2W + v 2 ), S v = ^a 2 p(v 2 ), C v = 0. (86) 
The kinematical backreaction vanishes, while the dynamical backreaction simplifies slightly to 

Vt> = + (2<S>d l d t V - (5,*) (d l -f + <9 l $)) . (87) 

The solutions to the linear system in Fourier space are given [33| by 

1 1 / x 2 

$ k = const., 5 = --$ k fcV, v = -$ k fc?7, a = a ( — 
6 3 V^o 

Noting as a result of (|69p that 

((a a $)(a Q $)) = - ($a a a Q $) , (89) 

one can find using these solutions that the effective energy density, pressure and equation of state of the backreaction 
from second-order perturbations are 

8ttG 2 _ r + 8ttG 2 _ 1 1 , nn ^ 

— «Pc ff = — 0, — ap cff = -Q, ^-^—^ (90) 



where 



r = |, (91) 
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and 

L=-j Wk) + 2$(k)) e - k x^ + 6 y^ (k)e - lk .x_|_ ; e~^ 2 |#)|VWy (92) 

are the linear and quadratic contributions respectively. (Recall that (5 = + (l/2)<5( 2 ), so the "linear" contributions 
include terms that are at second order in perturbation theory.) Here R = J W(x)y/hod 3 x/ J W(x)\/hd 3 'x. is the ratio 
of the spatial volume as evaluated on the background and the perturbed manifold, r then encodes the magnitude 
of straight averages across the perturbation compared to averages of the perturbations squared. In other words, it 
encodes the error implicit in the approach of, for example, Papers I and II, where it was assumed that L = 0. 7 It is 
not the aim of this study to evaluate the magnitude of r; however, some general statements may be made about the 
backreaction nonetheless. 

The effective pressure is positive definite. To accelerate the universe then requires a negative energy density, 
implying that r < —19. The backreaction then acts as a dark energy for w e g > —1/3, which requires r > —20. 
Therefore, if 

re (-20,-19) (93) 

then the backreaction arising in a pure dust universe will act to accelerate the spatial volume - as a dark energy, albeit 
one with a negative effective energy density which serves to reduce the Hubble rate. This is a restrictive requirement 
and argues against a backreaction that acts as a dark energy arising naturally in a pure dust universe. 

If instead, as seems likely, the average of linear modes vanishes on large scales due to their Gaussian nature, and 
the average of second-order modes is of the same order of magnitude as that from quadratic terms, we have 

0(r) M 1. (94) 

For r £ [—10, 10], we find that the effective equation of state lies between 

h < WeS < h (95) 

and in the case when the averages of second-order perturbations can themselves be entirely neglected we recover 

Wefl = (96) 

as found in Paper II. From this we can conclude that, for any system where 0(r) ps 1, the equation of state of 
backreaction from perturbations, up to and including second-order modes, is similar to warm dust and approximately 

WeS ~ — • (97) 
60 y ' 

Only for situations where r ~ —39/2 does the effective equation of state become such that the backreaction acts as 
a dark energy. This could be possible, for example, on smaller scales when the averages over linear modes cannot be 
neglected. 



IV. DISCUSSION 



In this paper we have formulated the averaged cosmological equations in the most general case for which the 
3+1 separation holds and considered the nonlinear regime of perturbation theory in an unfixed and in two simple 
gauges. The formalism we have derived is entirely general for any spacetime in any coordinate system, and with any 
combination of fluids. This formalism differs from that recently presented in [2311 as we have defined the Hubble rate 
from the change of a volume defined on an arbitrary surface, while the author of 23] defined his Hubble rate from the 
expansion scalar. Before specialising to a perturbative system, we also highlighted some issues with the interpretation 



7 While the second version of this paper was being prepared, the authors of |35|I studied the backreaction in Newtonian gauge including 
explicitly the averages of second-order perturbations and so effectively found r explicitly, in the case where ensemble averages of first-order 
perturbations vanish. 
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of the "backreaction" , depending on the choice of underlying system and of the lapse. These issues are illustrated with 
reference to perturbative models: in such systems the backreaction terms will in general include terms that otherwise 
appear in the background Friedmann equations; and for systems in a conformal time the energy densities should all 
be scaled by the conformal factor. It is therefore sensible to suggest that general statements concerning the "effective 
fluid" should be made with care, and that this effective fluid is best defined on a model-specific basis. 

One of the motivations for the development of this formalism, other than future study into fully general systems, 
was to consider the gauge issue in backreaction and, in particular, whether a choice can be found that renders the 
correction terms in their simplest form. To this end, we applied the general formalism to a Robertson- Walker universe 
perturbed to second-order in metric and fluid perturbations, working in an unspecified gauge. The resulting averaged 
equations can then be put into a required gauge by removing the redundant degrees of freedom. We then argued that 
the most obvious gauge choices are those that either remove or maximally simplify one or more of the backreaction 
quantities, with the simplest choices being the uniform curvature gauge, the conformal synchronous gauge, and the 
conformal Newtonian gauge. The system in uniform curvature gauge remains somewhat unwieldy despite the vanishing 
spatial curvature, with the kinematical and dynamical backreactions complicated instead with contributions from the 
shift vector. However, flat gauge does significantly simplify the use of Fourier modes and ensemble averaging. The 
system appears to reduce to a simpler form in the conformal Newtonian gauge, in which the anisotropic terms and 
the shift vanish and in which analytic solutions for simple cases are easily found. Working in conformal Newtonian 
gauge, we generalised the analytical estimates of Paper II to the case where contributions of the form ($) cannot be 
neglected, finding that, on large scales where the average of linear modes vanishes and that of second-order modes 
is of the same order as quadratic terms, the backreaction has an equation of state similar to warm dust. On smaller 
scales where the average of linear modes no longer vanishes it remains possible to argue for an equation of state more 
drastically different from dust, but the linear contributions would need to be in a narrow window more than an order 
of magnitude greater than those from quadratic terms. 
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